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ABSTRACT. 

We examine quantum properties of mesoscopic, Josephson coupled superconducting 
dots, in the limit that charging effects and quantization of energy levels within the dots 
are negligible, but quasi-particle transmission into the weak link is not. We demonstrate 
that quasi-particle resonances lead to current-phase relations, which deviate markedly 
from those of weak links connecting macroscopic superconductors. Results for the steady 
state dc Josephson current of two coupled dots are presented. 
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Recent theories of mesoscopic Josephson weak links [1,2] have revealed a variety of 
new quantization phenomena, which are absent from their macroscopic counterparts. 
For such systems, the size of the junction is less than or of order the phase breaking 
length, while the superconducting leads are macroscopic. In this Letter, we consider 
the quantum properties of weak links, formed when the superconductors themselves 
are mesoscopic. Linked mesoscopic, superconducting dots (LM dots), for which quasi- 
particles as well as Cooper pairs maintain phase coherence within the device, have been 
grown experimentally [3, 4], although no theory of their quantum properties currently 
exists. In this Letter, we demonstrate that the current-phase relation of LM dots is 
sensitive to resonant tunneling of quasi-particles into the junction from external, current 
carrying leads. This leads to a non-equilibrium distribution of quasi-particles within 
the weak link, which is missing from conventional descriptions of Josephson junctions[5]. 
Such resonances, which lead us to view LM dots, as electronic analogues of Fabry- Perot 
interferometers, carry a current of order the critical current and therefore significantly 
modify the transport properties of such structures. 

In a forthcoming publication [6] we shall present a complete theory capable of 
describing the steady state properties and slow dynamics of an arbitary number of LM 
dots. In this Letter we illustrate how scattering theory can be used to extract the key 
physics, by examining two such dots in one dimension, described by the Bogoliubov - de 
Gennes equation 

H(r)£(x) = , (1) 

where 

f[-(h 2 /2m)d 2 x + u(x)-ij] A(x) 
\ A*(x) -[-(h 2 /2m)d 2 + u(x) - fx] 

In this equation, fx is the condensate chemical potential, u(x) the normal scattering 
potential and A(x) the superconducting order parameter. A simple model of a pair a 
LM dots, shown in figure 1, is obtained by allowing A(x) and u(x) to be non-zero only 
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in regions of size L\ and L 2 , where the order parameter phase takes values 4>i an d 4>2 
respectively. Figure 1 shows a pair of LM dots connected by perfect, normal leads, to 
external reservoirs at chemical potentials (i a and fib and distinguishes the present approach 
from other descriptions [1,2, 5], where the sources of charge are of no consequence and the 
lengths L\ and L 2 are taken to be infinite. For LM dots, the system size L\ + L + L 2 
is assumed to be smaller than the quasi-particle phase breaking length and therefore a 
description, which incorporates quasi-particle phase coherence throughout the device is 
appropriate. 

To obtain such a description, consider the most general eigenstate of H belonging to 
eigen-energy E. In the regions where A(x) and u{x) vanish, this has the form 



/ Aexp[ikx] + L?exp[— ikx] 
I Cexp[iqx] + Dexp[— iqx] 



}&(x) = < 



A exp[ikx] + B exp[— ikx] 
C exp[iqx] + D exp[— iqx] 

A'exp[ikx] + B'expj— ikx] \ 
C'ex.-p[iqx] + D'exp[— iqx] J 



for x < -(Li +L/2); 
for-L/2 < x < L/2; 
for(L 2 + L/2) < x; 



(3) 



where h 2 k 2 /2m — \x = fi — h 2 q 2 /2m = E. In the absence of inelastic scattering, the 
quantum properties of such a structure can be described in terms of either a transfer 
matrix T or scattering matrix S, defined by 
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Once T is known, S can be constructed and vice versa[7,8]. Both T and S are functionals 
of all physical potentials, as well as functions of E. Since H is Hermitian, quasi-particle 
probability (though not charge) is conserved, and therefore S is unitary. 
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To describe the Josephson effect for such a structure, one must compute the current 
in the region — L/2 < x < L/2, as well as in the external leads. To this end we form the 
expectation value of the current density operator, with respect to the density matrix 
corresponding to incident distributions f"(E) of quasi-particles of type a, energy E 
along lead i, where a = +1(— 1) for particles (holes). In the simplest case, f^{E) 
would be a Fermi distribution, although the analysis outlined here can equally well 
describe transport properties arising from non-equilibrium (eg. hot electron) reservoir 
distributions. Expressions for the currents I r in the left and right leads respectively 
are written down in reference [9]. To obtain the current within the weak link, we introduce 
separate transfer matrices T 1: T 2 connecting plane wave amplitudes at x = to those 
in the left and right leads respectively and satisfying T = T 2 Ti. The wavefunction \I/(0) 
due to, for example, an incident particle from the left, is obtained by setting A = 1, 
D = B' = C = in equation (3) and using equation (4) to obtain the outgoing 
amplitudes B and C. The corresponding coefficients A", B" ', C" , D" , which we 
denote Ai + , C±+, -Di+, respectively are then obtained by acting on the vector 

(5b) with T\. In this way, by combining the action of S with T\ and plane wave 
amplitudes A^ a , Bi jCe , Ci^ a , Di jCe within the weak link, due to a quasi-particle of 
type a, incident along lead i can be constructed. For simplicity we restrict the present 
analysis to zero temperature, where the expectation value of the current density operator 



inside the weak link is found to be Ii n = I s + I qp , with 
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and 



I s = (2e/h) V r^{[|A-| 2 - \a-\ 2 ] + [q(E)/k(E)][\D i+ \ 2 - \C t+ \ 2 ]} 
i=i Jo 

I qp = (2e/h) ^dE{\A 1+ \ 2 -\B 1+ \ 2 + [q(E)/k(E)][\C 1+ \ 2 -\D 1+ \ 2 ]} 
Jo 

+ (2e/h) / dE{[k(E)/q(E)][\A 2 .\ 2 - |S 2 _| 2 ] + \C 2 .\ 2 - \D 2 _\ 2 } 
Jo 



(6) 



(7). 



Note that the right hand side of equation (6) represents the contribution from all occupied 
negative energy states, which has been transformed, using the particle-hole symmetry 
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relations A ia (E) = D*_ a (-E) and B ice (E) = —C*_ a (—E), to an integral over positive 
energies. The division of the total current into a sum of two currents is somewhat 
arbitrary. However if the reservoir potentials are equal, I qp vanishes, whereas I s may 
remain finite. Therefore in what follows, we refer to I s and I qp as the supercurrent and 
quasi-particle current respectively. 

Before proceeding, it worth noting that in deriving equations (6) and (7), the 
occupancy of incoming states from external reservoirs has been chosen such that the 
ground state expectation value of any local operator is preserved by the transformation 
from a closed to an open system. More precisely, consider a scatter of size V = L\ +L+L 2 
connected to leads of size L", which join together to form a closed system of size V + L" ', 
with periodic boundary conditions. If p is the density matrix, then the expectation value 
Tr pO(x) of an operator 0{x) can be evaluated using any convenient set of basis states. 
For a closed system the obvious choice is the set of eigenstates of H, satisfying periodic 
boundary conditions. However in the limit L" — > oo, such a choice is no longer useful 
and a trace over all incoming scattering states is preferred. It is important to note that 
if the incoming scattering states are interpreted as arising from external reservoirs, then 
certain properties are imposed on the reservoirs by the scatterer. For example at zero 
temperature, in the absence of a potential difference, expectation values are preserved 
only if all incoming quasi-particle states are populated; both incoming electron and 
incoming hole states. Such occupancies are non-intuitive, since a more natural choice 
of zero temperature reservoirs would perhaps populate only incoming electron states. 
The relative merits of different choices of incoming distributions will be discussed in 
a forthcoming publication[6]. Here we merely note that the existence of quasi-particle 
resonances is independent of such choice, although the detailed form of the current-phase 
relation does depend on the distribution of incoming states. 

As emphasized in reference [9] , due to non-conservation of quasi-particle charge, the 
chemical potential fi enters expressions for and I r explicitly and in steady state, must 

5 



be determined self-consistently by insisting that the currents are equal. In the present 
context one also notes that for an arbitrary phase difference = 4>i ~ 4>ii the internal 
current I m will not equal the current in the leads. Hence to obtain a dc Josephson effect, 
in which a current flows between superconducting dots of equal potentials, both fx and <f> 
must be determined from the steady state condition I in = I\ = I r . The solution to these 
equations yields a phase-current expression </>(/), which may be inverted to yield a more 
familiar current-phase relation I ((/)). If the current supplied by the external reservoirs is 
greater than a certain value, it may happen that no solution to these equations exists, 
in which case a critical current I c has been exceeded. Since the primary aim of this 
Letter is to highlight the role of quasi-particle resonances, full implementation of this 
self-consistent scheme will be relegated to reference[6]. One notes however that for the 
symmetric structures considered below, the equation I\ = I r is trivially satisfied with the 
choice [i = (Ha + A*b)/2. 

Except for a set of resonant energies, the quasi-particle contribution to L in is expected 
to be negligibly small, because for reservoir potentials less than typical values of |A(x)|, 
quasi-particle states decay on the scale of the superconducting coherence length, which 
in practice may be much greater than L\ and L2. For a clean junction, in the limit 
L l5 L 2 — > 00, these resonances correspond to bound state energies of the weak link 
considered by Bardeen and Johnson[10]. In the latter description, the penetration of 
quasi-particles into the weak link is ignored and therefore in the clean limit, the analysis 
presented here reduces to that of reference[10] in the absence of quasi-particle transmission 
through the superconductors. 

When solving equation (1), it is convenient to introduce a characteristic wavevector 
kp through the relation Ti 2 k 2 F /2m = \x and to divide both sides by [x. The resulting 
equation depends only on the dimensionless quantities E = E/li, u(x) = u(x)/lj,, 
A(x) = A(x)/(i and x = xkp. With this choice of scaling, equation (1) and therefore all 
scattering properties, do not depend explicitly on \x and it is convenient to express results 
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in terms of the dimensionless potential differences fi a = (fi a — fib = (/-* ~~ / /x. 

To illustrate the effect of resonances on the current-phase relationship of LM dots, figure 
2 shows plots of iin(0) versus for different choices of /it. Choosing |A(x)| = A in 
the superconducting regions of dimensionless length Li,L 2 , figure 2 shows results for 
the case L = 100, L\ = L 2 = 150, A = 0.04, u{x) = and an applied potential of 
P>ab = (a*o — A*6)/a* = 0.01. These values are typical of those attainable experimentally, 
where A is always much smaller than unity and the size of a dot is much greater than 
the coherence length £ = Aq 1 . The dotted, dashed and solid lines show results for 
fib = 0, flab/ 4 and /iab/2 respectively. The main part of the figure shows the internal 
current I\n{4>) = I s (4>) + Iqp{4>)-> while the insert shows the quasi-particle current. Except 
near the endpoints, the supercurrent I s increases almost linearly with phase, as <f> increases 
from —tv to 7r, reflecting the fact that under certain conditions, a phase difference across 
such a junction acts like a Galilean transformation[10]. In contrast I qp exhibits strong 
resonances over well defined intervals of <fi. In the limit L\,L,2 — > oo, when L » Aq 1 , 
the bound state energies of such a junction occur at [10] E n ~ ±<fi/L + (2n + 1)tv/L. As 
noted earlier these bound state energies become quasi-particle transmission resonances 
in a LM dots. For </> = ±7r, such a resonance occurs at E = and therefore quasi- 
particles from external reservoirs can enter the weak link. To understand these results 
in more detail, consider first the quasi-particle current shown by the dashed line, arising 
when fib = f-ab/^ and therefore fi a = 3/i a fe/4. As cj) increases from = — tt, the lowest 
resonant energy E Q first exceeds the highest incident quasi-particle energy fib from the 
right reservoir, leading to a decrease in the quasi-particle current by a factor of 1/2. 
As increase further, eventually the highest energy fi a of a quasi-particle from the left 
reservoir is exceeded and I qp switches off. For the case fi a = fib shown by the solid line, 
the quasi-particle currents from the two reservoirs switch off at the same phase, while for 
the case fib = 0, only the left reservoir contributes to I qp . 

It should be noted that the steps in 7 qp are of order the critical current through 
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the device, despite the fact that the superconductors are several coherence lengths long. 
This arises because resonant states within the weak link are formed from superpositions 
of particles and Andreev reflected holes, for which the currents add constructively. In 
a forthcoming publication, it will be shown that adding normal potential scatterers to 
the contacts between the external leads and superconductors does not significantly affect 
Andreev scattering within the weak link and therefore has only a marginal effect on the 
quasi-particle current, whereas potential scattering within the weak link suppresses both 
the supercurrent and the quasi-particle current. 

Figure 2 illustrates the effect of quasi-particle resonances for one value of the external 
potential difference only. For larger values of /i a — /it, more than one resonance can 
contribute and as shown by the solid lines of figure 3, the resulting internal current exhibits 
a non-trivial phase dependence. For completeness, the dashed lines show the external 
current in the leads and from the fact that these lines cross at several distinct values of 
4>, one concludes that the fully self-consistent solution will possess many branches. 

In this Letter, we have shown how current phase relations for LM dots can be obtained 
from a knowledge of the transfer matrices Ti, T<i and the associated scattering matrix 
S. The results obtained highlight the role of quasi-particle resonances in determining 
junction properties. Figures 2 and 3 demonstrate that the quantum properties of LM 
dots can differ markedly from those of more conventional Josephson junctions, a feature 
which should manifest itself in a range of junction properties. A key property of the 
current-phase relations obtained in this Letter is that Ii n (4>) is no longer an odd function 
of (f) and therefore its the phase average does not vanish. While one might be tempted to 
identify this phase average with a quasi-particle leakage current, it should be emphasized 
that it cannot be replaced by a phenomenological Ohmic term in a RSJ equation. Indeed 
in situations where the phase difference varies slowly with time, the quasi-particle current 
is transmitted as a series of pulses, reflecting the resonant nature of such electronic 
interferometers. 
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Figure Captions. 



Figure 1. 

The order parameter of a typical pair of LM dots, connected to external reservoirs 
at chemical potentials fi a and /ij,, with \x a > fi^. 
Figure 2. 

For fi a + fib = 0.01, this figure shows current-phase relations for a pair of LM dots of 
size Li = L 2 = 150, separation L = 100 order parameter A = 0.04 and u(x) = 0. The 
main body of the figure shows results for the total internal current ii n (</>) = I S (<P) + -^p (</>), 
arising when Jn, = (dotted line), /U a &/4 (dashed line) and /U a b/2 (solid line). The insert 
shows corresponding results for the quasi-particle current I qp (4>) only. 
Figure 3. For the same system used in figure 2, this figure shows the internal current 
-^in(</ ) ) (solid curves) and the external current / e xt(</ ) ) = h = Ir (dashed curves) for applied 
potential differences of magnitude ji a b = 0.02 (top left), 0.04 (top right), 0.06 (lower left) 
and 0.08 (lower right). 
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